Next-to-leading order QCD corrections to a heavy resonance 
production and decay into top quark pair at the LHC 

Jun Gao, Chong Sheng LiJ*| Bo Hua Li, and Hua Xing Zhu 

Department of Physics and State Key Laboratory of Nuclear Physics and Technology, 

Peking University, Beijing 100871, China 

C.-P-YuanS 

Department of Physics and Astronomy, 
Michigan State University, East Lansing, 48824, USA 
(Dated: August 16, 2010) 

Abstract 

We present a complete next-to-leading order (NLO) QCD calculation to a heavy resonance 
production and decay into a top quark pair at the LHC, where the resonance could be either a 
Randall-Sundrum (RS) Kaluza-Klein (KK) graviton G or an extra gauge boson Z' . The complete 
NLO QCD corrections can enhance the total cross sections by about 80% — 100% and 20% — 40% 
for the G and the Z', respectively, depending on the resonance mass. We also explore in detail 
the NLO corrections to the polar angle distributions of the top quark, and our results show that 
the shapes of the NLO distributions can be different from the leading order (LO) ones for the KK 
graviton. Moreover, we study the NLO corrections to the spin correlations of the top quark pair 
production via the above process, and find that the corrections are small. 
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I. INTRODUCTION 



The top quark is the heaviest particle so far discovered, with a mass close to the elec- 
troweak symmetry breaking scale, and closely related to various new physics models beyond 
the standard model (SM). Thus it provides an effective probe for the electroweak symmetry 
breaking mechanism and the new physics beyond the SM through studying its production 
and decay at colliders. The Large Hadron Collider (LHC) is running now with a center 
of mass energy yfs = 7 TeV, and will collect 1 fb _1 experimental data during the initial 
run. After this initial state the LHC will turn to y'i = 14 TeV, with a design luminosity of 
~ 10 fb _1 /yr there will be 8 x 10 6 top quark pairs and 3 x 10 6 single top quarks produced 
yearly. As a result of all these, the precision measurement of the top quark properties, 
such as the mass, the production cross sections, the kinematic distributions, and the spin 
correlation effects, will be one of the prime tasks in the experiments at the LHC, and any 
deviations from the SM predictions will definitely be a hint for new physics beyond the SM. 

To explore the connections between the new physics and the top quark, one possibility 
is to study the top quark pair invariant mass distribution and look for possible resonances 
since many new physics models predict the existence of a new resonance with a mass around 
TeV, which can decay into a top quark pair, such as the Technicolor [1], Topcolor {2], Little 
Higgs [3|], general Z' models [J, |5|], and Randall- Sundrum (RS) models [6|. In addition, 
in many such models the interaction between the heavy resonance and the top quark is 
enhanced as compared to the other fermions and the resonance will mainly decay into a top 
quark pair, for example, the Kaluza-Klein (KK) excitations of the graviton the gluon 8] 



as well as the weak gauge bosons 



9] in the extended RS models. Once we have discovered 



such a resonance in the top quark pair invariant mass distribution, the next step is to measure 



its spin and couplings, anc 
have been studied in Refs. 



finally determine the underlying new 
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phy sics dynamics, which 



111 ] . It has been suggested in Refs. [10|, ]ll| that it is possible 



to extract the spin and coupling information of the resonance from the top quark polar angle 
distributions and the spin correlations of the top quark pair. Those studies were carried out 
at the leading-order (LO) in QCD interactions. However, the next-to-leading order (NLO) 
QCD corrections may be large, for example, the NLO QCD corrections can enhance the 



cross sections of the single RS KK graviton or the Z' production by about 70% 
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13| 



and 20% 14j respectively, so it is necessary to examine whether the QCD corrections will 
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change the tree- level results and some conclusions of Refs. [10|, LL3J or not. In this paper 
we investigate the NLO QCD effects to a heavy resonance production and decay into a top 
quark pair, i.e., pp — » X(color singlet) — > ti, at the LHC, where the SU(3)c color singlet 
state X could be either a RS KK graviton G or an extra gauge boson Z'. 

The arrangement of this paper is as follows. Section [II] is a brief review to the relevant 
models. In Sec. II I II we show the details of the NLO calculations. Section [IV] contains the 
numerical results, and Sec. [V] is a brief summary. The Appendix collects some analytic 
results at the LO. 



II. THE MODELS 



A. The RS KK graviton 



In the RS model, a single extra dimension is compactified on a S 1 /Z2 orbifold with a 
radius r, which is not too large as compared to the Planck length. Two 3-branes, the Planck 
brane and the TeV brane, are located at the orbifold fixed points <fi = 0, 7T, respectively, and 
the spacetime between the two 3-branes is simply a slice of a five-dimensional anti-de Sitter 
geometry. The five-dimensional warped metric is given by 

ds 2 = e- 2kr ^ri^dx»dx u - r 2 d<\) 2 , (1) 

where <fi is the five- dimensional coordinate, and k ~ Mp is the curvature scale. By requiring 
kr ~ 12, one can suppress the Planck scale to Mpe~ hnr ~ O(TeV) on the TeV brane, and 
then solve the gauge hierarchy problem. The gravity fields are treated as fluctuations under 
the background metric, and after expanding the gravity fields in the extra dimension we get 



infinite massive KK gravitons, which can interact with the SM fields 15]. 








FIG. 1: Tree-level Feynman diagrams for the heavy resonance production and decay into a top 
quark pair. 



3 



The RS KK graviton can be produced through both the gg fusion and the qq annihilation 
at the LO as shown in Fig. HJ The detailed Feynman rules of the graviton couplings can be 
found in Re ? 



is given by 



12 



16| , and the propagator of the graviton in the unitary gauge in n dimensions 



with 
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where mx and Tx are the mass and the width of the heavy resonance, respectively. 



B. The extra gauge boson Z' 

1 — I 

The extra gauge boson Z' could arise from an additional U(l)' gauge symmetry [4|. It 
could also be the KK excitation of the electroweak gauge bosons. It can only be produced 
through the qq annihilation at the LO, and its generic couplings to quarks are as follow 

* QQ ~ 7 M I a L— — + a R 2 j ■ ( 4 ) 

We considered four combinations of at and aR, i.e., the pure vector coupling, ai = aR = 1, 
the axial- vector coupling, a^ = —aR = —1, the right-handed coupling, ol = 0, or = 1, 
and the left-handed coupling, a L = 1, a R = 0, which are denoted by Z[, Z' 2 , Z' z and Z' 4l 
respectively. The propagator of the Z' in the unitary gauge is given by 

PZ ' = k 2 -m\ + imx Tx + W ' 

In our calculations of the process pp — > X (color singlet) — > ti, what we mainly concern 
about are the ratios of the NLO results to the LO ones, so it is not necessary to specify the 
actual values of all the couplings. We simply assume the mass of the heavy resonance to be 
around TeV scale, which is not yet excluded by the current experiments. Besides, we only 
consider the narrow resonance cases and fix Tx/rnx = 1% at both the LO and the NLO. 
We do not expect our conclusions to be largely modified even if Tx/rax increases to be at 
the order of 10%. Detailed discussions on the SM backgrounds and the discovery potential 



of the process can be found in Refs. 
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III. THE NLO FORMALISM 



The complete NLO QCD corrections to the process pp — » X (color singlet) — > ti can be 
factorized into two independent gauge invariant parts, i.e., the heavy resonance produced at 
the NLO with a subsequent decay at the LO, and produced at the LO with a subsequent 
decay at the NLO, similar to the cases studied in Ref. jr/J . The box diagrams, and the 
corresponding real correction diagrams, that connect the initial and the final states do not 
contribute to the squared matrix elements up to the NLO as the heavy resonance is a SU (3)c* 
color singlet particle. This whole procedure can be illustrated as follows 

\m^ 2 \ 2 = \M t ;::\ 2 ®\Mt e :?®\Px\\ 

M ^ M %* = ^\M^\ 2 9 (M^M^ + \^4^\ 3 9 (M^M^)} 9 \P X \ 2 , (6) 

we have suppressed the possible Lorentz indices here for simplicity. 

We calculate the full squared matrix elements using the propagators (Px) of the heavy 
resonance given in Sec. [Til which can incorporate the full spin correlations between the 
production and decay processes in order to generate the correct kinematic distributions. We 
carry out all the QCD calculations in the 't Hooft-Feynman gauge and use the dimension 
regularization scheme [3] (with the naive 75 prescription [jij) in n = 4 — 2e dimensions to 
regularize all the divergences. The one-loop Feynman diagrams for the production and the 
decay of the heavy resonance are shown in Fig. [2J 

Fig. |3] shows the real correction Feynman diagrams for the production and the decay of 
the heavy resonance. The infrared divergences of the real corrections are extracted by using 



the two cutoff phase space slicing method [201 ] . Due to the limited space here we do not 
reproduce the details of the method. 

IV. NUMERICAL RESULTS 

In our numerical calculations we choose the input parameters to be m top = 171 GeV, 



mz = 91.188 GeV, and a s (mz) = 0.118 21_|. The running QCD coupling constant is 



evaluated at the three-loop order J2l| with rif = 5, and the CTEQ6M (CTEQ6L1) parton 
distribution function (PDF) set [22] is used through the NLO (LO) calculations. We set 
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FIG. 2: Some one-loop Feynman diagrams for the production and the decay of the heavy resonances. 
Others not shown can be obtained by the exchange of the external quark or gluon lines. 




FIG. 3: Some real correction Feynman diagrams for the production and the decay of the heavy 
resonances. Others not shown can be obtained by the exchange of the external quark or gluon 
lines. 

both the renormalization and factorization scales equal to the mass of the heavy resonance, 
unless specified. Besides, in the two cutoff phase slicing method there are two arbitrary 
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cutoff parameters, i.e., the soft cutoff S s and the collinear cutoff 8 C . We have checked the 
cutoff dependence of all our numerical results, and found that the dependence is negligibly 
small for 5 S < 1 x 10~ 3 , so we choose 5 S = 1 x 10~ 3 and 5 C = 5 s /50 to obtain the numerical 
results presented below. 



A. The total cross sections 



In Fig. H]we show the NLO K factor, defined as the ratio of the NLO cross section <Jnlo 
to the LO cross section ctlo, as a function of the heavy resonance mass at the LHC with 
different center of mass energies. We can see that the total NLO QCD corrections can be 
large, which can enhance the total cross sections by about 80% — 100% and 20% — 40% for 
the G and all four types of Z' bosons, respectively, depending on the resonance mass. The 
NLO corrections from the production part are dominant, while the ones from the decay part 
are relatively small, but can still reach above ten percent in some regions. Our results of 

nn 

the NLO K factors of the production part agree with the ones given in Refs. [13|, UJ, where 
the total cross sections have been summed over the spins of the heavy resonance directly. 
In the following parts of our paper we will only show the results of the total NLO QCD 
corrections for simplicity. We further present the ratios of the total cross sections from the 
different channels for the graviton at both the LO and the NLO in Fig. It can be seen 
that the contribution from the gg channel is dominant at the low mx value region due to 
the large PDF of the gluon, and the contribution from the qq channel becomes important 
at the high mx value region since the PDF of the valence quark decreases more slowly than 
the gluon. And the NLO corrections can change the ratio of the contribution from the qq 
channel to the one from the gg channel significantly. 



B. The polar angle and invariant mass distributions 



It has been shown in Refs. 
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ll| that the polar angle distributions of the top quark 



are the key points to extract the spin and coupling information of the heavy resonance. 
The definition of this polar angle depends on the reference frame and axis chosen, here we 
considered two kinds of polar angles, one is the Collins-Soper angle 8s 23], which is defined 
to be the angle between the top quark momentum and the axis that bisects the angle between 



7 




1000 1500 2000 2500 3000 1000 1500 2000 2500 3000 



m x (GeV) m x (GeV) 



FIG. 4: The NLO K factors as functions of the heavy resonance mass at the LHC, the solid and 
dotted lines correspond to including the total QCD corrections and the QCD corrections from the 
production part alone, respectively. The four groups of the curves from the top to the bottom 
correspond to the G, Z' 2 , Z^(Z' 4 ), and Z[, respectively. 
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FIG. 5: The ratios of the total cross sections from different channels for the graviton as functions 
of the graviton mass at both the LO and the NLO. 

the momentums of the incoming hadrons (j>a and —ps) in the ti rest frame, and for the Z' 
we can define 9* as the angle in the ti rest frame between the top quark momentum and the 
incident quark momentum which can be determined by the longitudinal boost direction of 
the ti rest frame at the LHC [lo| . 

In Figs. [MS we show the normalized polar angle distributions of the top quark at the 



8 




FIG. 6: The normalized top quark polar angle distributions at the LO and the ratios of the 
normalized NLO distributions to the LO ones, at the LHC with ^/i = 14 TeV for mx = 800 GeV. 




FIG. 7: The normalized top quark polar angle distributions at the LO and the ratios of the 
normalized NLO distributions to the LO ones, at the LHC with y/s = 14 TeV for mx = 1500 GeV. 

LO and the ratios of the normalized NLO distributions to the LO ones at the LHC with 
yfs = 14 TeV. At the LO, we can use both the cos 9s and cos#* distributions to distinguish 
the Z' and the G as their distributions have significantly different shapes. At the same time 
we can also use the cos 8* distribution to distinguish between the Z[ 2 and the Z' z 4 since the 
latter ones have a large forward-backward asymmetry. Furthermore, the differences of the 
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FIG. 8: The normalized top quark polar angle distributions at the LO and the ratios of the 
normalized NLO distributions to the LO ones, at the LHC with ^/s = 14 TeV for mx = 3000 GeV. 

polar angle distributions between the Z[ and the Z' 2 are very small for mx around 1 TeV 
or heavier, thus it is not possible to separate them through the polar angle distributions. 
We present the LO squared helicity amplitudes in the Appendix, which can explain the 
behavior of the LO polar angle distributions. After including the NLO corrections, we can 
see that for all the Z ! the changes of the distributions are negligibly small, which are no 
more then a few percent. But for the G, as the increasing of the resonance mass the NLO 
corrections can change the shapes of the distributions, for example, the NLO corrections 
make the distributions decrease more quickly at the both ends and can reach about 10% 
for nix — 1500 GeV. The corrections can be as large as 30%, and change the shapes 
of the distributions significantly for mx = 3000 GeV, which do not change the fact that 
the distributions of the G and the Z' are greatly different. Note that for a heavy enough 
resonance, the top quarks produced are highly boosted, which means the decay products 
of the top quark are close to each other and form a top jet. Recently, several methods 



based on the jet substructures have been proposed [241 ]. which may be used to detect such 
a top jet efficiently, so it is possible to measure the NLO QCD effects to the polar angle 
distributions for a graviton with a mass of several TeV at the LHC. The reason that the NLO 
distributions for the graviton differ from the LO ones is that the NLO corrections change the 
ratio of the contributions from the gg and qq channels, as shown in Sec. IIV A] which have 
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different shapes of distributions. As the resonance mass increases, the contributions from 
these two channels become comparable, so the changes are more significant. We further 
study the scale and PDF uncertainties of the NLO polar angle distributions for a graviton 
with mx = 3000 GeV. As shown in Fig. [HI the uncertainty from the scale dependence is 
negligibly small, and the PDF uncertainty is within 10%, which is still small as compared 
to the NLO corrections. Here, we use two more PDF sets in the NLO calculations, i.e., the 



MRST2004nlo 



25| and MSTW2008nlo [26| PDF sets. 




FIG. 9: The scale and PDF uncertainties of the NLO polar angle distribution for the graviton at 
the LHC with y/s = U TeV and m x = 3000 GeV. 

In Fig. [10] we present the invariant mass distributions of the top quark pair including 
the NLO QCD corrections. At the LO they are just the Breit-Wigner distributions with a 
center value mx and a width Fx- While at the NLO the heavy resonance can decay into a 
top quark pair plus a hard gluon, so the NLO corrections increase the distributions in the 
lower invariant mass value region, and the changes of the distributions are more significant 
as the resonance mass increases. We also studied all the above distributions at the LHC 
with y/s = 7 TeV, and the results are similar. 



C. The spin correlations 

One of the unique features of the top quark is that it decays before the strong interaction 
can depolarize its spin. Thus, it is possible to extract the spin information of the produced 
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FIG. 10: The normalized top quark pair invariant mass distributions at both the LO and the NLO 
at the LHC with ^/i = 14 TeV. 

top quark by studying the angular distributions of the decay products. For a spin up top 
quark (or a spin down anti-top quark), the decay angular distribution of the zth decay 
product is given by [271 ] 

I dT 1 

- (1 + aicosxi) , (7) 



r T d(cosXi) 2 

where i could be quarks, b,u,c,d,s, or leptons, Ui,l; x% is t ne angle between the zth decay 
product and the spin quantization axis in the top rest frame, and at are the correlation 
coefficients. For the charged leptons, a\ = 1 exactly, which means the charged leptons are 
maximally correlated with the top spin direction. 

The spin correlations of the top quark pair also can be used for the identification of the 
heavy resonance, but the precision measurement of them is more difficult at the LHC. In 
order to study the spin correlations of the top quark pair production, the following double 
differential cross section is usually considered, 

~~Ti TT37 =T = 7 ( 1_ Aa i a j cos xt cos Xj + h a i cos xt + cos xj) , (8) 

a d(cosxT)d(cosXj ) ^ 

neglecting the interference between the top spins we have 

<r(t t t t + J - afoii + *jj t ) 



A 



o-(t t t t + J + cr{t^t i + t^) ' 
q(t t t t + tfa) - a(t^ + tjtj,) 

a(t^ + t^) - er(£ t t t + t^) 



<r(t t t t + t&) + a{t^ + tifr) ' 
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In our following calculations we use the helicity basis in the tt center of mass frame, which 
means xtixj) is defined to be the angle between the t(t) direction in the tt center of mass 
frame and the corresponding decay product direction in the t(t) rest frame. 




1000 1500 2000 2500 3000 



m x (GeV) 

FIG. 11: The top quark pair spin correlation coefficients at the LO as functions of the heavy 
resonance mass at the LHC with yfs = 14 TeV. 

In Fig. [11] we show the LO results of the top quark pair spin correlation coefficients 
as functions of the heavy resonance mass at the LHC with yfs = 14 TeV. According to 
symmetry analysis we have 

A(Z' 3 ) = A(Zi), b±(G,Z[,Z' 2 ) = 0, 

b + (z' 3 ) = b-(z' 3 ) = -b + (z' 4 ) = -b4z' 4 ), (io) 

and at the LO for the axial vector Z 2 , 

A(Z' 2 ) = -1, (11) 

which can be seen from the helicity amplitudes in the Appendix. With the increasing heavy 
resonance mass, all the coefficients in Fig. [TT] approach —1 due to the fact that the cross 
sections for the ti with the same helicities vanish as the heavy resonance mass goes infinity. 
In Table [TJ, we list some typical NLO results of those coefficients. We can see that the NLO 
QCD corrections are rather small, about 1% — 2%, and can be neglected at the LHC. We 
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mass(GeV) 


resonance 


A(LO) 


A(NLO) 


MLO) 


6+ (NLO) 


MLO) 


6_(NLO) 




G 


-0.783 


-0.778 


















-0.832 


-0.826 














800 


z' 2 


-1.000 


-0.986 


















-0.904 


-0.896 


0.947 


0.943 


0.947 


0.943 




Z'a 


-0.904 


-0.896 


-0.947 


-0.943 


-0.947 


-0.943 




G 


-0.933 


-0.913 
















z[ 


-0.949 


-0.933 














1500 


z' 2 


-1.000 


-0.980 


















-0.974 


-0.955 


0.986 


0.977 


0.986 


0.977 




Z'i 


-0.974 


-0.955 


-0.986 


-0.977 


-0.986 


-0.977 



TABLE I: The top quark pair spin correlation results at the LHC with y/s = 14 TeV. 

also investigate the cases for yfs = 7 TeV, and the results are almost the same at both the 
LO and the NLO. 

V. CONCLUSIONS 

We have calculated the complete NLO QCD corrections to a heavy resonance production 
and decay into a top quark pair at the LHC, where the resonance could be either a RS KK 
graviton G or an extra gauge boson Z' . Our results show that the total NLO K factors can 
reach about 1.8 — 2.0 and 1.2 — 1.4 for the G and all four types of Z' bosons, respectively, 
depending on the resonance mass. And the NLO corrections from the production part are 
dominant, while the ones from the decay part are relatively small but can still reach above 
ten percent in some parameter regions. We also explore in detail the NLO corrections to the 
polar angle distributions of the top quark, and our results show that the NLO distributions 
are almost the same as the LO ones for all four types of Z' bosons, while the shapes of the 
NLO distributions can be significantly different from the LO ones for the G, depending on 
the mass of the resonance. Moreover, the NLO corrections can also change the shapes of the 
top quark pair invariant mass distributions. Finally, we study the NLO corrections to the 
spin correlations of the top quark pair, and find that the corrections are negligibly small. 
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Appendix 



\M\\ 



\M\ 



\M\< 



Asm 4 (9/2)A, 
4cos 4 (6/2)A, 

4/3 2 sin 4 (fl/2)„4, 
4/3 2 cos 4 (#/2)„4, 



(12) 



(13) 



(14) 



In this appendix we give the individual nonvanishing LO squared helicity amplitudes for 
a heavy resonance production and decay into a top quark pair, qq (gg) — > X — > ti . For all 
the Z' mediated processes, 

' (1 - f3 2 ) sin 2 (6)A, for helicities {+ - ++}, {+ }, 

{- + ++} and {- + } 

for helicities {H h} and { — h H — } 

for helicities {H 1 — } and { — I h}, 

for helicities {H h} and { — h H — } 

for helicities {H 1 — } and { — I h}, 

(1 - (3 2 ) sin 2 (#)„4/4, for helicities {+ - ++} and {+ } 

(1 - (5f sin 4 (^/2)^l, for helicities {+ - -+} 
(1 + /3) 2 cos 4 (^/2)^, for helicities {+ - +-}, 

(1 - f3 2 ) sin 2 (^)^/4, for helicities {- + ++} and {- + } 

(1 - f3) 2 sin 4 (^/2)^, for helicities {- + +-} 
(1 + /3) 2 cos 4 (^/2)^, for helicities {- + -+}, 

and for the graviton mediated processes through the qq annihilation and the gg fusion, 

/3 2 (1 - f3 2 ) sin 2 (2#)£/64, for helicities {+ - ++}, {+ }, 

{- + ++} and {- + } 

/3 2 (1 + 2 cos(#)) 2 sin 4 (#/2)£/16, for helicities {+ - -+}, 

and {- + +-} 
f3 2 (l - 2 cos(#)) 2 cos 4 (#/2)£/16, for helicities {+ - +-}, 

and {- + -+}, 

3/3 2 (l - f3 2 ) sin 4 (#)£/128, for helicities {+ - ++}, {+ - 

{- + ++} and {- + } 

3/3 2 sin 6 (#/2)cos 2 (#/2)£/8, for helicities {+ - -+}, and {- + +-} 
3f3 2 sin 2 (#/2) cos 6 (fl/2)£/8, for helicities {+ - +-}, and {- + -+}, 



\M\ 2 



(15) 



\M\ 



G,qq 



(16) 



\M\ 2 



G,gg 



(17) 



with 



A 



(s - m|,) 2 + m|,r|,' 
16 



B 



s — m 2 G ) 2 + m 2 G Y 2 G '' 



(18) 



where 9 is the polar angle between the momenta of the top quark and the light quark (or 
gluon) in the center of the mass frame of the top quark pair, s is the square of the center 



of mass energy, /3 = yl — Amf op /s, and the squared amplitudes have been summed and 
averaged over the color of the external particles. 
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